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Abstract With the increasing interplay between ex¬ 
perimental and computational approaches at multiple 
length scales, new research directions are emerging in 
materials science and computational mechanics. Such 
cooperative interactions find many applications in the 
development, characterization and design of complex 
material systems. This manuscript provides a broad 
and comprehensive overview of recent trends where pre¬ 
dictive modeling capabilities are developed in conjunc¬ 
tion with experiments and advanced characterization 
to gain a greater insight into structure-properties re¬ 
lationships and study various physical phenomena and 
mechanisms. The focus of this review is on the intersec¬ 
tions of multiscale materials experiments and modeling 
relevant to the materials mechanics community. After a 
general discussion on the perspective from various com¬ 
munities, the article focuses on the latest experimental 
and theoretical opportunities. Emphasis is given to the 
role of experiments in multiscale models, including in¬ 
sights into how computations can be used as discovery 
tools for materials engineering, rather than to “simply” 
support experimental work. This is illustrated by exam¬ 
ples from several application areas on structural mate¬ 
rials. This manuscript ends with a discussion on some 
problems and open scientific questions that are being 
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explored in order to advance this relatively new field of 
research. 
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Multiscale modeling • Integrated computational 
materials Engineering (ICME) • 3-D microstructure 
characterization • Characterization methods 

1 Introduction 

Recent advances in theoretical and numerical meth¬ 
ods, coupled with an increase of available high perfor¬ 
mance computing resources, have led to the develop¬ 
ment of large-scale simulations in both Materials Sci¬ 
ence and Mechanical Engineering, with the goal of bet¬ 
ter characterization and optimization of materials per¬ 
formance. These advanced computational capabilities 
extend over a large span of length scales—ranging from 
the atomistic to the continuum scale—and serve as in¬ 
vestigative tools to get a greater insight into structure- 
properties relationships. As these predictive modeling 
capabilities become more comprehensive and quanti¬ 
tative, a comparable level of details from experiments 
and characterization tools is now not only necessary 
to validate these multiscale models, but also to mo¬ 
tivate further theoretical and computational advances. 
The present complexity of the multiscale approach leads 
also to the need for new algorithms to bridge length 
scales and for model reduction. Strong coupling be¬ 
tween experimental data and numerical simulation is 
required for predicting macroscopic behavior with a 
fine description of local fields (i.e. mechanical, compo¬ 
sitional, electrical, etc.). An illustration of such syner¬ 
gies is the coupling of experimental tools such as Trans¬ 
mission Electron Microscopy (TEM) and synchrotron- 
based X-ray microscopy that enable the collection of 
microstructural statistical data comparable to the out- 
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put of three-dimensional (3-D) Crystal Plasticity Finite 
Element Method (CPFEM)-based simulations [1,2]. 

As illustrated in Fig. 1, based on a citation report 
extracted from ISI Web of Knowledge, there is an in¬ 
creasing trend of articles in the fields of both Materi¬ 
als Science and Mechanics to report on both multiscale 
modeling and experiments. In Fig. la, we consider those 
that denote a topic area of Coupling Experiments and 
Modeling^f A similar trend is observed for articles con¬ 
taining Multiscale Modeling^' in the title for research 
areas limited to the same fields of Materials Science 
and Mechanics. Both progressions in the number of 
journal articles published highlight the increasing cross¬ 
pollinations between the modeling community and the 
experimental community. While the term “ coupling^' is 
far too inclusive to extract any real correlation between 
the emergence of advanced characterization methodolo¬ 
gies and growth of research into ‘'^Coupling Experiments 
and Modeling^', it is interesting to put in perspective 
this increase of publications with respect to the succes¬ 
sive introduction of novel experimental characterization 
methodologies such as Scanning Electron Microscopy 
(SEM) [3], Atom Probe Tomography (APT) [4], SEM- 
based Electron Back Scattered Diffraction (EBSD) [5] 
and 3-D high energy synchrotron X-ray approaches [6]. 
As such, the increasing use of various individual exper¬ 
imental techniques within modeling paradigms is illus¬ 
trated in Fig. lb. 

1.1 Background 

Just what is this emerging field? After all, many mul¬ 
tiscale modeling strategies discussed in the mechan¬ 
ics and Materials Science communities couple modeling 
and experiments intrinsically, i.e., they involve mod¬ 
els of different physics at multiple scales and exper¬ 
imentally validate those models at the corresponding 
length/time scales. But because of the interdisciplinary 
nature of the fields, a novel branch in Materials Sci¬ 
ence that relies on the synergy between experimen¬ 
tal and computational approaches at multiple length- 
scales is opening up new frontiers and research direc¬ 
tions at the crossroads of both traditional Computa¬ 
tional and Experimental Materials Science, in addition 
to the emerging field of Integrated Computational Ma¬ 
terials Engineering (ICME). Such cooperative interac¬ 
tions find many applications in the development, the 
characterization and the design of complex material 
systems. Yet more possibilities are waiting to be ex¬ 
plored and many new questions of physical, numerical, 
analytical or characterization nature have materialized. 
These things all make this emerging field an extremely 
fruitful and promising area of research. 


Almost all problems in Materials Science are multi¬ 
scale in nature: fundamental small scale physical mech¬ 
anisms and processes occurring at the atomic length 
scale and femtosecond time scale have a profound im¬ 
pact on how materials perform at larger spatial and 
time scales. Despite this inherent multiscale charac¬ 
ter, effective macroscopic models and measurements are 
used—often with satisfactory accuracy—that account 
intrinsically for the effects of the underlying microscopic 
processes. For example, a conventional tensile test di¬ 
rectly measures the ultimate tensile strength, the maxi¬ 
mum elongation, the reduction in area, and determines 
a phenomenological uniaxial stress-strain curve. Such 
an experiment can be used effectively to calibrate a 
Hollomon-type power law relation between the stress 
and the amount of plastic strain regardless of the fun¬ 
damental microstructural mechanisms responsible for 
the macroscopic response. 

However, the above class of macroscopic models and 
experiments have some limitations. The first obvious 
limitation is the complete neglect of microscopic mech¬ 
anisms that are sometimes of interest. In our tensile 
test example, we may want to know how changes in mi¬ 
crostructure would lead to a change in the mechanical 
response. It is often of interest to know and charac¬ 
terize the microstructural defects such as dislocations, 
martensite formation, and mechanical twinning that are 
at the origin of plastic deformation, not just the macro¬ 
scopic flow. Such information requires modeling and ex¬ 
perimental tools that can track such defects. These in¬ 
clude, for example, CPFEM or discrete dislocation dy¬ 
namics (DDD) on the modeling side and TEM, X-ray 
tomography or SEM-based EBSD on the experimental 
side. Another limitation is associated with the heuris¬ 
tic nature of the macroscopic models and experimen¬ 
tal measurements. These effective models often have an 
empirical foundation limiting their predictive capabili¬ 
ties, while the interpretation of experimental data may 
be challenging in the case of complex material systems 
subjected to the intricate interplays of thermodynam¬ 
ics, transport (diffusion and phase transformations) and 
non-linear behavior. As we will describe later, one of 
the advantages of detailed multiscale models is their 
ability to switch which physical mechanisms are and 
are not active in the material to decouple and exam¬ 
ine the importance of different mechanisms and help 
in the interpretation of experimental data. Finally, an¬ 
other limitation is that it is difficult to gain insight into 
the margins and uncertainty due not only to test-to- 
test variation, but also the natural stochasticity at the 
lower length scales of engineering materials that makes 
an “exact” knowledge of the underlying microstructure 
and active mechanisms impossible. 
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On the topic of "Coupling Experiments and Modeling" 
■MB With "MttlHscale Modeling" in the title 



Year 



(a) Papers with “ Coupling Experiments and Modeling" (b) Papers focused on modeling and various experimen- 

as topic area and with “Multiseale Modeling" in the title. tal techniques. 

The circle symbols denote the introduction of specific 
experimental technologies. 


Fig. 1: Evolution of the number of journal articles published in the fields of Materials Science and Mechanics 
according to ISI Web of Knowledge. 


For these reasons, one might be tempted to switch 
completely to a detailed microscopic characterization 
(either model-wise, experimentally or both) that has 
the highest resolution and strongest physical founda¬ 
tion. However, this is not always a favorable strategy; 
not only because the microscopic models and experi¬ 
mental setup are often difficult to handle, but also be¬ 
cause it requires complex procedures to extract the in¬ 
formation of interest and perform meaningful data min¬ 
ing. 

This is where the coupling between experiments and 
simulations is essential. By combining multiscale mod¬ 
eling and corresponding experiments, one hopes to take 
advantage of the simplicity, efficiency and mechanistic 
insight gained from the models, as well as the phys- 
icality, meaningful reproducibility and “reality check” 
provided by the experiments. 

1.2 Different communities, different perspectives 

Various perspectives from diverse communities lead to 
different combinations of experiments and simulations. 

Perspective of Solid Mechanics: As pointed out 
by Hortstemeyer [7], the Solids Mechanics community 
tends to adapt a “top-down” or ^^downscaling" approach. 
Such tendency is driven by the need to characterize ma¬ 


terials performance through the use of internal state 
variables (ISVs) and degrees of freedom (DoF) repre¬ 
senting an underlying microstructure at the pertinent 
length scales. Stemming for the most part from the Mi¬ 
cromechanics community [8-10], this led to the develop¬ 
ment of hierarchical and nested methods in which both 
numerical/theoretical models and experiments probe 
properties and mechanisms down to various subscales 
to characterize the evolution of ISVs and thus describe 
the materials behavior at the macroscale. For exam¬ 
ple, a large body of work is currently being conducted 
on the quantification and characterization of the be¬ 
havior of small size systems (i.e. when the structure 
size compares with the microstructure size) such as 
micro- and nano- electromechanical systems (MEMS 
and NEMS) [11-14] or 3-D printing structures with thin 
walls and beams [15-17]. 

Perspective of Materials Science: The Materials 
Science community focuses generally on a “ bottom-up" 
or ^‘upscaling" approach. Indeed, the basis of Materi¬ 
als Science involves studying the structure of materials 
and their associated length scales, and relating them 
to their properties. While Materials Science researchers 
(notably those in the Mechanics community, who study 
the interplay between microstructure and mechanical 
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response) have always connected theory and experi¬ 
ments, it was not until the mid-nineties and the emer¬ 
gence of advanced characterization tools (see Fig. 1) 
that multiscale studies have proliferated to investigate 
structure-property relations. Synergies between exper¬ 
iments and modeling emerging from this community 
are motivated by the rationale that the different struc¬ 
tures composing a material (e.g. grain size and texture, 
dislocation network, inclusions) dictate its performance 
properties. For example, many studies now focus on 
the characterization of the microstructural evolutions 
through in situ observations [18,19] or on 3-D statisti¬ 
cal descriptions of microstructure [20-22]. 

Perspective of Applied Mathematics: The above- 
mentioned communities require the proven methodolo¬ 
gies of the Applied Mathematics community [23, 24] 
supporting the characterization and study of materi¬ 
als behavior across multiple length scales both from 
the modeling and experimental aspects. The major fo¬ 
cal points of this community lies in the scale bridging. 
Thus, methodologies such as dimensional analysis [25, 
26], fractal analysis [27-29], statistical analysis [30,31], 
graph theory [32,33], parallel computing [34-36], and 
model reduction [37,38] started from the applied math¬ 
ematics community and were successfully used by both 
modelers and experimentalists in the Materials Science 
and Solid Mechanics communities. 

Although the combination of experiments and simu¬ 
lations is common to many disciplines, this review will 
focus on recent progress made in metals. The authors 
readily recognize many areas where the integration of 
experiments with multiscale modeling has affected the 
Materials Science community, for example in the area 
of polymers [39,40], or biomaterials [41-43]. The pur¬ 
pose of this manuscript is to provide the Materials Sci¬ 
ence community with a coherent overview of the status 
of recent progress made into how Computational Ma¬ 
terials Science can be exploited as discovery tools for 
materials engineering and directly integrated with ex¬ 
perimental work. It is our hope that this summary will 
help the reader to understand (i) the recent successes 
of predicting various physical phenomena and mecha¬ 
nisms in materials systems enabled by the collaboration 
between experimentalists and modelers; (ii) the main 
obstacles that one has to face when experiments are in¬ 
tegrated into various computational mechanics frame¬ 
works; and (iii) a perspective on the upcoming issues 
to solve in order to make this emerging field of Mate¬ 
rials Science a more powerful and predictive approach. 
The manuscript is organized as follows. Section 2 high¬ 
lights recent synergies at various scales ranging from 
the atomistic scale to the continuum scale from an ex¬ 


perimental perspective. Section 3 provides a discussion 
on the recent developments of models and upscaling 
strategies between the various length scales. Section 4 
discusses the roles of experiments in multiscale models. 
The review ends with a discussion on some problems 
that have to be addressed in order to successfully lever¬ 
age the experiment-modeling synergy. 

2 Recent evolutions in experimental techniques 
across multiple length scales 

The quantitative description of microstructure is an 
essential requirement for multiscale models describing 
the macroscopic properties of materials [44]. Thus, the 
rapid evolution of experimental characterization tools 
opens up new opportunities to measure a wide variety 
of local fields accurately within a large volume and sub¬ 
sequently incorporate these measurements into a model 
at the corresponding scale. Not only can this type of ex¬ 
perimental data be used to calibrate models, it can also 
drive the development and improvement of new theo¬ 
retical approaches. 

Recent developments of experimental instrumenta¬ 
tion including (i) new capabilities for mechanical mea¬ 
surements and microstructural analyses at the nanome¬ 
ter length scale (Section 2.1); (ii) the development of 
full field measurements for different parameters beyond 
displacement field quantification (Section 2.2); and (iii) 
the improvement of 3-D characterizations (e.g. tomog¬ 
raphy; Section 2.3), have a major impact on the mod¬ 
eling community. 

2.1 Nanoscale microstructural and mechanical analysis 

The advent of optimized materials systems engineered 
at the near-atomic scale requires the quantification of 
microstructural features (local composition, crystalline 
structure) and local mechanical measurements at that 
same scale. 

Local chemical composition measurements: Over 
the years, a large number of small-scale measurement 
techniques have been developed to quantify the local 
chemical composition in structures through recent ad¬ 
vances in analytical High-Resolution (HR)TEM [45,46] 
and in analytical Scanning (S)TEM. For instance, new 
STEM systems employ very high sensitivity X-ray de¬ 
tectors [47, 48] allowing for very sensitive 2-D and 3- 
D measurement of composition. Coupled with electron 
energy-loss spectroscopy and high-angle annular dark- 
field imaging, it is now possible to reach a spatial accu¬ 
racy of few atoms, allowing for the local identification of 
the chemical composition and, in specific cases, the ox¬ 
idation state [49]. However, the acquisition of 3-D char¬ 
acterization of the composition within a TEM remains 
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challenging. For this, the Materials Science community 
turns to APT [49]. 

Even though it is exercised by a small community, 
APT [4,50“52] has shown the most progress in improved 
resolution of small length scale compositional measure¬ 
ments. The technique enables the identification of the 
chemical nature of atoms using time-of-flight mass spec¬ 
trometry and, at the same time, their positions within 
the sample with sub-nanometer resolution. Data ob¬ 
tained from APT provide new insights about mecha¬ 
nisms controlling the evolution of local chemical com¬ 
position. Figure 2 exemplifies a recent study of grain 
boundary segregation in an Al-Mg alloy [53] using APT 
to reveal the mechanisms controlling dynamic precipita¬ 
tion and segregation during severe plastic deformation 
of aluminum alloys. Due to the relatively small dimen¬ 
sions 10® nm®) of the materials volumes and the 
reasonable number of atoms 10^ atoms) that can be 
characterized using this technique, an increasing num¬ 
ber of studies couple atomistic simulations and APT 
(ISI Web of Knowledge citation report extracted 460 
articles between 2002 and 2015, see Fig. lb). 

In situ mechanical field measurements: Experi¬ 
mental techniques being used to measure the elastic 
fields at the nanometer length scale are now also avail¬ 
able. Among them, high resolution X-Ray Diffraction 
(XRD) [54-57] is particularly interesting since it allows 
for lattice orientation, and lattice strain tensor coor¬ 
dinates to be measured with great accuracy. Because 
the full lattice strain tensors are available, correspond¬ 
ing mean stress tensors may be calculated directly from 
the experimental data at the sub-grain level unambigu¬ 
ously using classical elasticity theories. Such measure¬ 
ments are especially useful to study strain/stress evolu¬ 
tion in situ and conduct residual stress analysis. They 
are well-suited to validate polycrystalline models such 
as CPFEM [58-60] or to inform local mechanistic crite¬ 
ria such as crack emission or damage at grain boundary, 
for example. Thus, Van Petegem et al. [58] conducted in 
situ mechanical testing using XRD on nanocrystalline 
Ni to observe the evolution of inter-granular stresses 
and used molecular dynamics (MD) simulations and 
crystal plasticity models to help with the interpreta¬ 
tion of the experimental results. Similarly, recent fem¬ 
tosecond XRD experiments have been combined with 
atomistic simulations to characterize lattice dynamics 
in shock-compressed materials to study the ultrafast 
evolution of microstructure [61,62]. 

Mechanical properties measurements at small 
scales: Motivated by the development of nanoscale sys¬ 
tems (e.g. NEMS) and by the explosion of sub-grain 
level simulations (21,684 publications reported from ISI 






Fig. 2: Atom probe analysis of Mg concentration at a 
grain boundary (GB) triple junction in a Al-Mg alloy 
deformed at 200°C: (a) Al atoms; (b) Mg map with 
indication of the two concentration profiles shown in 
(c) for Profile A, and (d) for Profile B. Depletion and 
segregation of Mg is quantified close to the GBs within 
a 2.5 nm step size (the sampling volume thickness is 
1 nm). Reprinted from Ref. [53], Gopyright (2014), with 
permission from Elsevier. 


Web of Knowledge for the 2010-2015 time period in the 
research area solely limited to the field of Materials Sci¬ 
ence on the topic of ^‘‘Molecular Dynamics'\ for exam¬ 
ple), mechanical testing at the micro-scale is currently 
a very popular field of research. 
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A conventional approach to obtain local mechani¬ 
cal properties is through hardness measurements. Some 
development in modeling the nanoindentation process 
is given in Section 3.2, but sufficed to say that this 
popular experimental technique (8,952 publications re¬ 
ported from ISI Web of Knowledge for the 2010-2015 
time period in the research area solely limited to Mate¬ 
rials Science) enables the extraction of various quanti¬ 
tative characteristics of the materials being tested from 
the load-displacement curve recorded. These include, 
for example, the elastic behavior (both during loading 
and unloading), macroscopic yield stress, hardness, lo¬ 
cal contact characteristics, and internal stresses in near¬ 
surface layers [63]. The magnitude of the indentation 
force applied and the shape of the indenter change not 
only the characteristic size of a locally deformed region 
but also the contributions of the elastic and plastic de¬ 
formations and in the sub-surface active deformation 
mechanisms. Therefore, when coupled with such exper¬ 
iments, various types of computational models are se¬ 
lected depending on the relative deformation mecha¬ 
nisms activated: 

— At the initial stages of the indenter penetration, 
when the nucleation and migration of individual 
structural defects come into action, atomistic sim¬ 
ulations are used to mimic nanoindentation exper¬ 
iments. These support the analysis and interpreta¬ 
tion of the experimental results to gain insight on 
the deformation mechanisms [64-66]. For example. 
Lodes et al. [64] studied the indentation size effect 
and pop-in behavior in CaF 2 single crystals, using 
both nanoindentation experiments and MD simula¬ 
tions. 

— For later stages of indenter penetration when con¬ 
siderable plastic deformation occurs, atomistic tools 
are no longer appropriate, rather CPFEM-type cal¬ 
culations become more relevant simulation tools [67, 
68]. For example, by combining microindentation 
with a CPFEM approach, Yao et al. were able 
to quantify the plastic parameter of slip systems in 
tungsten single crystal [69]. For polycrystalline ma¬ 
terials, local behavior measured by nanoindentation 
was combined with a corresponding CPFEM model 
to study the anisotropic response of a-titanium [67]. 

In situ nano-hardness experiments within a TEM [70- 
73] consisting of the integration of nanoindentation with 
real-time electron imaging also merit mentioning. Clas¬ 
sically, the main objective of these experiments is real¬ 
time study of the interplay between deformation mech¬ 
anisms and the mechanical response at the nanome¬ 
ter length scale. Such experiments are computationally 
well-suited for atomistic modeling simulations [74, 75] 
or discrete dislocation dynamics simulations [76-78], 


where these computational paradigms are used to pro¬ 
vide new insights on the fundamental physical phenom¬ 
ena dominating dominating the mechanical response. 
For example, Bufford et al. [78] used in situ nanoin¬ 
dentation to study plasticity and work hardening in alu¬ 
minium with incoherent twin boundaries. Discrete dis¬ 
location dynamics simulations were used to clarify the 
role of stacking-fault energy in the difference of defor¬ 
mation mechanisms between nanotwinned copper and 
twinned aluminium. 

Novel experimental advances at the nanoscale allow 
for a precise insight of the compositional, microstruc- 
tural and mechanical states of materials. As indicated 
by multiple examples above, this type of characteriza¬ 
tion is of primary importance to obtain quantitative 
information for the analysis of active mechanisms (e.g. 
the plastic behavior and defect activities near a grain 
boundary). However, in order to obtain accurate and re¬ 
liable measurements at this length scale, numerical sim¬ 
ulations have to be consistent with this scale to identify 
relevant mechanisms, or obtain materials properties by 
reverse modeling. 

2.2 Local field measurements and analysis 

Recent experimental advances in local field measure¬ 
ments is another area of Materials Science where exper¬ 
iments have been impactful within the modeling com¬ 
munity since such experiments have enabled a higher 
degree of verification and validation of models beyond 
the macroscopic and average behavior. 

Surface displacement fields from microgrid tech¬ 
niques: Although not very accurate and with a low 
resolution, photoelasticimetry and holographic interfer¬ 
ometry are simple and useful legacy tools to visual¬ 
ize the displacement fields at surfaces. For example, 
in 1987, Boehler and Raclin [79] leveraged the natu¬ 
ral roughness of a tensile sample to develop a technique 
of pseudo relief in order to deduce from two subsequent 
images a complete displacement field on the surface. 
Since then, refinements in the resolution and visualiza¬ 
tion of the surface displacement fields came from print¬ 
ing a regular microgrid on tensile SEM samples for in 
situ experiments [80-82]. Initially used for deep drawing 
tests, millimeter circles were printed on the surface of 
the sample before deformation to (i) measure the direc¬ 
tion of the principal deformation axes and the principal 
strains and (ii) visualize the occurrence of the necking 
in a complex stamped part. The microgrids were usu¬ 
ally made out of two orthogonal sets of parallel lines 
of sputtered metal onto the sample surface. The thick¬ 
ness of the lines was a few tenths of a micrometer and 
mesh size was a few square micrometers. One of the 
advantages of using regular grids like this is easy de- 
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tection of shear localizations [83]. More recently such 
a technique has been coupled with CPFEM simula¬ 
tions [84-86] in order to understand the correlations 
between the mechanical behavior and the underlying 
microstructure. CPFEM calculations are carried out on 
a mesh of the microstructure based on orientation imag¬ 
ing microscopy (generally EBSD mapping) and the sim¬ 
ulation results are compared to the experimental local 
strain field measurements obtained from digital image 
correlation of the microgrid. Few attempts have been 
made so far to get similar in-volume information (see 
Section 2.3). 

Local field measurements through Digital Image 
Correlation: The development of displacement field 
measurements has been expanding with the automa¬ 
tion of quantitative image analysis techniques [87-92] 
such as Digital Image Correlation (DIC). DIC allows to 
go beyond the use of a regular grid and avoids some bias 
due to microgrids regularity and the fraction of surface 
covered by the grid lines. Speckle techniques are now 
regularly used to improve the accuracy of strain mea¬ 
surements at the grain level and at finer scales where 
plastic strain localization is manifested, for example, 
as physical slip bands [2,92]. Stinville et al. [92] stud¬ 
ied damage accumulation processes in nickel-based su¬ 
peralloy using a DIC approach with speckle patterns 
achieved by etching. By combining DIC and EBSD (see 
Fig. 3) they were able to detect low levels of strain and 
strain heterogeneities at the microstructural features. 

The popularization of EBSD measurements: As 

seen in the examples above, microgrid and DIC tech¬ 
niques are generally combined with EBSD measure¬ 
ments. The development of EBSD in the last decade has 
allowed displacement field measurements to be comple¬ 
mented with crystallographic orientation field measure¬ 
ments [93]. The automation processes and the reduction 
of probe size have enabled the mapping of fairly large 
surface area containing many grains with a spatial res¬ 
olution now in the order of 0.1 /rm [94]. EBSD measure¬ 
ments recorded during the deformation lead not only to 
the knowledge of the crystallographic texture evolution, 
but also to information about heterogeneities [88, 94- 
98]. As described later on in Section 3.1, such informa¬ 
tion is used to reconstruct “synthetic” equivalent repre¬ 
sentative microstructures. Average misorientation met¬ 
rics, such as the Kernel Average Misorientation (repre¬ 
senting the numerical misorientation average of a given 
pixel with its six neighbors within an EBSD image), 
have been developed to quantify the correlation be¬ 
tween plastic deformation and EBSD data. Orienta¬ 
tion gradients and geometrically necessary dislocation 
(GND) densities can be calculated directly from the ex- 



Fig. 3: Tensile (x direction represents the loading di¬ 
rection) strain field e^x (l^ft), and shear strain field Cxy 
(center) from DIC with images at high magnification 
after 0.5%, 1%, 1.5%, 2% and 3.5% macroscopic strain 
in tension. Associated EBSD map (before loading) and 
SEM images at the different macroscopic strains are re¬ 
ported on the right. Reprinted from Ref. [92], Copyright 
(2015), with permission from Elsevier. 


perimental measurements [99-102] and compared with 
non-local crystal plasticity formulations [103,104], for 
example. Additionally, although not as popular yet, sta¬ 
tistical analysis [105-107] (85 publications reported by 
ISI Web of Knowledge for the 2010-2015 time period in 
the research area solely limited to Materials Science on 
the topic of “EBSD” and “statistical analysis”) on large 
EBSD datasets are being used to extract and study 
quantitative relationships between material microstruc¬ 
tural features and the mechanical properties. One of the 
main limitations with EBSD measurements lies in the 
degradation of the signal for highly deformed materials. 
Pixel indexing becomes difficult and the accuracy of the 
measurement is limited. Nevertheless, it is possible to 
overcome this difficulty by using the pattern quality in¬ 
dex, which is another way to investigate the microstruc¬ 
ture. Such a technique has been applied, for example, 
to recrystallization problems [108]: recrystallized grains 
are easily indexable while the non-recrystallized do¬ 
mains appear with a very low quality index. 
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Local field measurements through TEM-based 
approaches: For very fine microstructures and/or lo¬ 
cal measurements, Automated Crystallographic Orien¬ 
tation Mapping obtained via TEM (ACOM-TEM) is 
an efficient tool [109, 110]. Well-adapted for severely 
deformed samples (by equal channel angular extrusion, 
for instance) and for nanograin-size materials, this tech¬ 
nique, similar to EBSD, allows for the collection of spe¬ 
cific information on the mechanisms at the grain level 
and to confirm the model prediction at the same scale. 
Finer spatial resolution of local misorientation down to 
a few nanometers is now accessible through the use of 
Precession Electron Diffraction (PED-ACOM). For ex¬ 
ample, Taheri et al. [Ill] accurately estimated dislo¬ 
cation density in a pure oxygen-free high thermal con¬ 
ductivity copper deformed by rolling and annealed at 
different temperatures from PED data using the Nye 
tensor. In addition to the popularity of these automated 
crystallographic orientation mapping techniques, other 
TEM-based approaches also enable accurate measure¬ 
ments of local fields a very fine scales. For example, 
HRTEM measurements of atomic plane displacements 
have been used to explore ferroelectric-ferroelastic in¬ 
teractions in PbTiOs [112]. The development of dark- 
field electron holography is another example of a tech¬ 
nique used to map the strain field at the nanoscale [113]. 
A recent publication [114] reviews the limitations and 
advantages of various experimental methods for mea¬ 
suring strain at this scale. The emergence of such tech¬ 
niques opens up new opportunities within the modeling 
community to validate and verify models. 

2.3 3-D measurements and in-volume characterization 

Surface observations, even during in situ experiments, 
have proven inadequate to quantify microstructure evo¬ 
lutions and mechanical field measurements, especially 
for subsurface behavior. Planar sections through a vol¬ 
ume can reveal information about the microstructural 
state, but local mechanical information can only be de¬ 
duced through numerical simulations. Several options 
are available as alternatives to the displacement field 
measurements techniques described in the previous sec¬ 
tion. The advent of such techniques proved valuable 
within the modeling community. 

3-D reconstruction: Serial sectioning, often associ¬ 
ated with EBSD measurements, allows for the recon¬ 
struction of the crystalline structure of a Representative 
Volume Element (RVE) and the associated grain orien¬ 
tations. This type of technique is widely used in con¬ 
junction with CPFEM simulations [115-117]. Recently, 
the fastidious tasks of repetitively polishing, drift cor¬ 
rection and SEM observations have been replaced by 
in situ focused-ion beam (FIB) layer removal combined 


with EBSD measurements [2,118-120]. For example, 
microstructure reconstructions permit a way to evalu¬ 
ate, in the case of multiphase materials, the phase con¬ 
nectivity and associated influence on the macroscopic 
behavior and damage/fracture processes (Fig. 4) [117]. 


Fig. 4: 3-D phase distribution in a Mo-TiC composite. 
The top micrographs represent X-ray tomography of 
the TiC phase distribution (a), and the Mo-TiC com¬ 
posite (b). The specimen size is 30x30x30 /rm^. The 
bottom figures are composite reconstruction of 3-D se¬ 
quential EBSD maps using a FIB for successive sec¬ 
tioning: (a) TiC phase, (b) Mo-TiC composite. The an¬ 
alyzed volume is 30x30x6.6 ^m^.Reprinted from [117], 
Copyright (2012), with permission from Elsevier. 


While the experimental approach is well established, 
the 3-D reconstruction and the fine-scale definition of 
grain and phase interfaces are still difficult. In fact, as 
the phase and crystal orientation is recorded pixel by 
pixel, the resulting interfaces are irregular and com¬ 
posed of small fiat rectangular facets [95]. This is par¬ 
ticularly detrimental for the accuracy of numerical cal¬ 
culations for phenomena addressing grain and/or inter¬ 
phase boundary effects on the stress-displacement field. 
The regular orientation of the facets with respect to 
the main macroscopic directions and the existence of 
sharp edges impact strongly the local estimation of the 
mechanical fields. This is of importance for simulation 
of any process in which interfaces play a major role 
(e.g. intergranular fatigue damage, grain boundary dif¬ 
fusion, nanostructured materials). Efforts are presently 
dedicated to smooth the surfaces obtained through 3-D 
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reconstruction techniques with numerical approaches, 
such as the level set technique, to overcome defects of 
both experimental digitalization and discrete numerical 
meshing [121]. 

In-volume characterization: XRD and X-ray com¬ 
puted tomography provide allow for in-volume charac¬ 
terization [122-131]. The idea is fairly simple: X-rays 
are used to penetrate bulk metallic samples and the lo¬ 
cal stress/strain state is deduced from the evolution of 
the peak location and its enlargement—revealing the 
evolution of the microstructure. A virtual representa¬ 
tion of the microstructure is constructed and modeled 
using a model (most commonly CPFEM) to simulate 
the evolution of the internal state variable mediating 
the polycrystal. Simulations are then compared directly 
to experimental diffraction data. For example, Odd- 
ershede and coworkers [124] measured the stress field 
around a notch in a coarse grained material using XRD 
and compared their measurements with a continuum 
elastic-plastic finite element model. The measured and 
simulated stress contours were shown to be in good 
agreement except for high applied loads, correspond¬ 
ing to experimentally observed stress relaxation at the 
notch tip. As illustrated in Fig. 5, lattice strain pole 
figures (SPF) combined with CPFEM-based methodol¬ 
ogy [128] are also used for quantifying residual stress 
fields within processed polycrystalline components. 

In the past. X-ray tomography has been limited to 
thin samples with a low density of defects, however with 
the development of near-held High-Energy Diffraction 
Microscopy (nf-HEDM) and the increasing availability 
of synchrotron sources, measurements of volumes on 
the order of a few cm^ are now possible. For exam¬ 
ple, this technique has recently been used to experimen¬ 
tally track the evolution of internal damage [56] or as a 
means to provide microstructural input (grain structure 
and crack morphology) for a CPFEM model embedded 
in a coarse homogeneous FE model [132]. As pointed 
out by Whithers and Preuss [125], such experimental 
techniques, especially when combined with mesoscale 
models (e.g. CPFEM), are becoming increasingly pop¬ 
ular for characterizing local phenomena such as fatigue 
crack growth or damage evolution. X-ray tomographic 
techniques are also used to obtain a complete 3-D char¬ 
acterization of architectured materials and track their 
evolution during deformation [133]. Under certain con¬ 
ditions, in situ measurements can be made, allowing for 
direct correlation with the numerical simulation results. 
XRD experiments combined with modeling capabilities 
offer a promising means to quantify mechanical fields 
but are so far limited by the time required to set up 
and run such experiments. 
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Fig. 5: High energy X-rays combined with CPFEM: 
(Top) Lattice strain pole figures from experiment and 
simulation, (bottom) The difference between the experi¬ 
ment and simulation strain pole figures. Reprinted from 
Ref. [128], Copyright (2014), with permission from El¬ 
sevier. 


2.4 Toward 4-D characterization 

Recent advances now not only enable 3-D measure¬ 
ments (see Section 2.3) but they also integrate a fourth 
dimension—time—and track temporal evolution at res¬ 
olutions as brief as picoseconds. Such characterization 
has been coined “four dimensional (4-D) Materials Sci¬ 
ence.” Notably, techniques such as hybrid environmen¬ 
tal scanning electron microscopy (ESEM)-STEM to¬ 
mography [134,135] and dynamic TEM [136-138] are 
emerging as high time resolution electron microscopy 
techniques with different variants enabling the real-time 
visualization of microstructure evolution. It is antici¬ 
pated that atom-probe tomography will gain this capa¬ 
bility in the future [4,139-142]. A demonstration of the 
ultrafast temporal resolution achieved so far is shown 
for example in Fig. 6 [143]. The balance between space 
and time accuracies are well illustrated by in situ X-ray 
tomography evolutions [56]. This allows in situ record¬ 
ing of solidification, recrystallization, and phase trans¬ 
formation kinetics in polycrystalline materials. 

The new opportunities offered by the experimental 
techniques and examples described in Sections 2.1-2.3 
illustrate how synergies between experiments and mod¬ 
eling can facilitate microstructure analysis (all the way 
down to the nanoscale). Motivated by a constant need 
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Fig. 6: Dynamic TEM images of the time-evolving showing the growth dynamics of radial dendritic grains in 
germanium. Reprinted with permission from Ref. [143]. Copyright (2012), with permission from Elsevier. 


to characterize deformation mechanisms more precisely, 
some improvements and iterations on the accuracy of 
both the experimental methodologies and correspond¬ 
ing simulations are necessary to evaluate the effects 
of microstructural variation on material response crit¬ 
ically. The non-exhaustive list of 3-D and 4-D charac¬ 
terization techniques described above will require ma¬ 
jor scientific advances and improved resolutions (both 
spatial and temporal); from new modeling and exper¬ 
imental approaches to new instrumentation. However, 
as we will discuss in Section 5, given the volume of 
data that can be acquired, such progress will need to 
be accompanied by (big) data management in order to 
efficiently and effectively collect, format, analyze, store, 
mine, and correlate large data sets linking process, mi- 
crostructure, and materials properties [144]. 

3 Recent advances in modeling 

The previous section highlighted recent experimental 
developments and illustrated how such measurements 
were integrated within corresponding models. In this 
section, we review examples of modeling and associ¬ 
ated numerical techniques using available experimen¬ 
tal information across many length scales including (i) 
modeling addressing the effect of heterogeneous mate- 
rials(Section 3.1); (ii) the development of approaches 
to model nanostructured devices and hybrid materials 
(Section 3.2); and (iii) the emergence of experimental 
measurement techniques incorporated into numerical 
frameworks directly (Section 3.3). 

3.1 Modeling of heterogeneous materials 

Understanding links between micro- and macro- scale 
properties is the key to predicting the behavior of media 
with a heterogeneous microstructure. As illustrated in 


Heterogeneous medium Representative 




Equivalent homogeneous medium 
Fig. 7: Mean-field versus full-field approach. 


Fig. 7, two general approaches using microstructural 
information exist in the literature: 

Mean-field approaches, initiated by Eshelby [8], are 
based on the solution of the inclusion of an ellipsoid 
(representing the grain) in the equivalent medium (rep¬ 
resenting the heterogeneous material, such as a poly- 
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crystal). They have been improved by, for example, 
taking into account elastoviscoplastic strains [145], or 
the influence of various microstructure parameters such 
as the presence of different phases [146] and the grain 
size distribution [147]. The main advantage of such ap¬ 
proaches is their computational cost-effectiveness allow¬ 
ing microstructures of more than 1,000 crystalline ori¬ 
entations to be randomly generated and tested numer¬ 
ically; however, since the strain is assumed constant 
within each grain, this technique cannot account for in- 
tragranular heterogeneities and fine descriptions of the 
microstructure, such as real grain shapes and orienta¬ 
tions. 

Full-field approaches have been developed to better 
take into account the latter specificities, by introduc¬ 
ing relevant RVEs [148]. The simulation costs, however, 
severely limit the size of the simulation cell. Moreover, 
the question of the size of the RVE and its numerical ho¬ 
mogenization remain open problems for polycrystalline 
materials: several multiscale strategies focusing on dif¬ 
ferent numerical schemes can be used, e.g., multilevel 
FEMs, for which Schroder [149] gives a recent review, 
or so-called “equation-free” approaches, for which the 
Heterogeneous Multiscale Method (HMM) [24] is well- 
known. 

Microstructural representations are required in ei¬ 
ther approach and are obtained either by using mi¬ 
crostructures obtained from image-based experimental 
characterization techniques (see Section 2.1 on using 
TEM images. Section 2.3 on using X-ray-based mea¬ 
surements, and Section 4.2 on using experiments as 
model input in general) or by generating a “synthetic” 
statistically representative microstructure numerically 
with varying degrees of randomness [150-152]. 

While appealing, exact numerical replicates of ac¬ 
tual microstructures in models are nonetheless cumber¬ 
some and necessitate appropriate algorithms to process 
such big data. Many models use relevant epresentative 
volume elements (RVEs) as an alternative to determine 
the global macroscopic behavior of a heterogeneous ma¬ 
terial. The general approach is to use a small set of 
microstructural descriptors (e.g. determined from 2-D 
experimental cross-sections) covering features includ¬ 
ing material composition, dispersion, and microstruc¬ 
ture morphology. From this, statistically equivalent 3- 
D microstructures are then generated [153-155]. As il¬ 
lustrated in Fig. 8, a good example is the recent mi¬ 
crostructure reconstruction algorithms based on a par¬ 
tial set of spatial correlations inspired by the field of 
solid texture synthesis [155]. Recent works stemming 
from the modeling and applied mathematics commu¬ 
nities are oriented towards numerical methods allow- 


2D Characterization 
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neighborhoods from each 2D 
scan. This is done for every 
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Fig. 8: A schematic description of overall process for 
microstructure reconstruction using partial set of de¬ 
scriptors derived from 2-D characterizations. Reprinted 
from Ref. [155], Copyright (2016), with permission from 
Elsevier. 


ing for the automated conversion of experimental data 
to entities usable by EE simulations. For example, an 
anisotropic tessellation technique [156] improved the 
reconstructions of actual microstructures from EBSD 
images, whereas an adaptive orientation reconstruction 
algorithm to get 3-D volumes [157] has been proposed 
for nf-HEDM images. Similarly, Lieberman et al. [158] 
suggested an improved technique for determining the 
normal vectors to grain boundaries extracted from nf- 
HEDM acquisitions, allowing a much finer prediction of 
void nucleation sites. 

It should be noted however that, while synthetic 
RVEs based on experimental statistical data have been 
considered historically (e.g. Ref. [159]), the limited size 
of such microstructures as being representative of the 
overall behavior may lead to outliers and often prove 
to be insufficient to obtain accurate local field informa¬ 
tion [129]. 

Homogenization and effective behavior of het¬ 
erogeneous materials: Many theoretical results ex¬ 
ist as far as the homogenization of random media con¬ 
cerned: starting with linear elliptic equations, these re¬ 
sults predict the existence of a constant deterministic 
tensor called “effective tensor”, that yields a displace¬ 
ment field close to the solution of the original prob¬ 
lem [160-162]. However, the actual computation of the 
value of this effective tensor remains an issue, since its 
approximation converges more slowly as the size of the 
domain is increased [163]. Some recent alternatives ex¬ 
ist in order to circumvent this drawback [164]. It is also 
possible to go further and address a stochastic equiv- 
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alent homogeneous model [165-167], even if this has 
not been specifically applied to polycrystalline materi¬ 
als modeling yet. 

There have been some attempts in the literature to 
truly couple probabilistic models, such as the multi¬ 
scale stochastic FEM by Xu et al. [168]. Other strate¬ 
gies are based on specific coupling schemes: whereas a 
surface coupling [169] was introduced to deal with lo¬ 
calized uncertainties, specific applications of a superpo¬ 
sition method, called the Arlequin method [170], have 
been proposed recently [171,172]. 

In addition to these multi spatial scale problems, 
it is also possible to address the simulation of physi¬ 
cal phenomena exhibiting different well-separated time 
scales, such as damage evolution for fatigue loading. 
For the latter case, the adaptation to time of classi¬ 
cal schemes of space homogenization is straightforward 
and leads to strong reductions in computational costs 
when compared with the initial full-scale problem: for 
example. Fuel and Aubry [173] give a tentative review 
of the so-called periodic time homogenization method. 
It is possible to deal with several scales simultaneously 
in time and space, e.g. Ref. [174]. Similar adaptations 
can be proposed in the stochastic framework: for exam¬ 
ple, HMM can be used in MD [175] or in kinetic Monte 
Carlo schemes [176]. 

3.2 Mechanical modeling of nanoscale systems 

Modeling and simulations strategies described in Sec¬ 
tion 3.1 are often challenged when the volume or the 
scale at which they are performed are of the same or¬ 
der as the material microstructures they must consider. 
Three domains where this is the case are considered 
here. 

Nanostructured materials: The low-dimensionality 
of nanostructures is constrained enough that continuum 
concepts are not applicable. Simulation tools such as 
atomistic modeling become necessary to simulate the 
elastic and mechanical properties [177]. Similar con¬ 
straints exist for nanosystems, such as NEMS [178,179] 
or nano-layered systems used for opto-electronics [180]. 
However coupling and comparison to corresponding ex¬ 
periments at this length-scale are challenging [181] due 
the general discrepancy in loading rates and time scales 
experimentally measurable and the ones achievable nu¬ 
merically. 

Architectured or hybrid materials: This class of 
materials motivates a lot of research in gaining ma¬ 
terials properties not achievable through “monolithic” 
materials [182,183]. As such, using hybrid structures 
to expand the property space results in the develop¬ 
ment of new numerical and experimental multi-physics 
approaches. However, there is a wide-range of technical 


challenges yet to be addressed. Among these is the prob¬ 
lem of similarity of length scale associated the struc¬ 
tural topology of such materials and the one associated 
with the physical phenomenon that such material is try¬ 
ing to address. This is exemplified by gradient materi¬ 
als [184] or by nanostructural devices fabricated by 3-D 
printing. In some cases, the part thickness is such that 
the material is locally a single- or a multi-crystal with 
its own anisotropy and properties variations from one 
location to another [185,186]. Similar concerns hold in 
nanostructured materials in which the typical size of 
the microstructure is comparable to the length scale of 
deformation mechanisms (such as the deformation and 
fracture of nanocrystalline materials [187] and heavily 
deformed pearlite wires [188]). Again, atomic calcula¬ 
tions are generally used in these systems. 

Local investigation of the microstructure are par¬ 
ticularly important for nanoindentation measurements 
(discussed in Section 2.1), where the atomic structure 
of the material has to be accounted for [77]. Thus, MD 
or DDD calculations are used to estimate the material 
behavior below the indenter [189] or during micro pil¬ 
lar compression [190]. Figure 9 shows a typical example 
of an experimental quantification of the displacement 
field in a copper micropillar; the measurement is de¬ 
duced from Die analysis of a speckle pattern deposited 
by FIB [190]. Combinations of EBSD and displacement 
maps allow the local properties to be deduced by in¬ 
verse calculation. Similar trends are noted for the anal¬ 
ysis of HRTEM observations. A complete simulation at 
the atom and dislocation level is needed to quantify 
experimental observations [191]. The main difficulties 
are then to define a representative volume that is large 
enough to include the local variation of the material, 
but small enough to limit the calculation size and time. 
Parallel calculations and model reduction are developed 
to enhance the simulation capacity. Nevertheless, the 
amount of deformation is still rather limited. The def¬ 
inition of the boundary conditions is very challenging 
and this problem is so complicated that a multiscale 
strategy is not usually possible at present. 

3.3 Simulated images as a complement to experimental 
characterization 

Well-established experimental methods such as electron 
and X-ray diffraction are being implemented into com¬ 
putational models as a complement for the interpreta¬ 
tion of experimental results. Simulated images can be 
done just like in real environments in which the model 
environment replicates the actual experimental settings 
and measurements acquisition. The goal of such mod¬ 
eling is two-fold: (i) verify the fidelity of the assump¬ 
tions and results obtained by a given model and (ii) 
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Fig. 9: Map of localized slips along Xi during com¬ 
pression of a copper micropillar. Reprinted from [190], 
Copyright (2014), with permission from Elsevier. 



Fig. 10: (Top) Diffraction contrast STEM of disloca¬ 
tions: imaging and simulations. Reprinted from [192], 
Copyright (2011), with permission from Elsevier. 
(Bottom) Comparison between virtual and exper¬ 
imental SAED patterns produced by (left) simu¬ 
lated and (right) experimentally observed in a-Al 203 
1711(1 0 1 1)(1 0 1 l)interface. Reprinted (adapted) 
from [193], Copyright (2015), with permission from El¬ 
sevier. 


supplement experimental studies to reveal convoluted 
microstructural information collected by a given exper¬ 
iment. 

Virtual diffraction [193,194] patterns generated by 
atomistic simulation, is a good example of this. These 
models use classical algorithms from kinematic diffrac¬ 
tion theory to create both selected-area electron diffrac¬ 
tion (SAED) patterns and XRD line profiles within an 
atomistic simulations. While this idea has been around 
for quite some time in Materials Science [195-197], the 
recent advances in virtual diffraction methods lie in the 
predictability of the methodology by performing virtual 
diffraction during an atomistic simulation via explicit 
evaluation of the structure factor equation without a 
priori knowledge of the grain boundary structure (in 
contrast to legacy studies that computed the structure 
factor equation from known structures). Recent exam¬ 
ples include studies of grain boundary structures [193] 
(see Fig. 10), microstructural evolution during deforma¬ 
tion [198,199], and deformation mechanisms [194,200]. 

Another great example of the replication of experi¬ 
ments within a computational framework is the STEM 
simulated image [192,201-203] (see Fig. 10). Image sim¬ 
ulations for the bright field and dark field TEM imag¬ 
ing uses a numerical algorithm based on the scattering 
matrix formalism within the context of the column ap¬ 
proximation. For each image pixel, the algorithm com¬ 
putes a full convergent beam diffraction pattern. Such 
simulated images have been shown to be useful for crys¬ 
talline defect analysis for a variety of diffraction config¬ 
urations. Note that image simulation of TEM disloca¬ 
tion images has a long history dating back to the 1960s, 
however the recent novelty lies in the application to the 
particular mode of scanning TEM imaging rather than 
conventional TEM and again also in the predictability 
of the method rather than a mere confirmation of the 
experimental characterization. 

4 The role of experiments in multiscale models 

Figure 11 highlights the information flow between ex¬ 
periments and models. The National Research Coun¬ 
cil recognized three distinct purposes for experiments 
in ICME efforts: (i) classic validation showing whether 
the models are consistent with what can be measured; 
(ii) model input that can set, e.g., model scaling pa¬ 
rameters or microstructure; and (iii) collection of large 
amounts of data that can be mined [205]. This section 
explores each of these in turn, but McDowell and Ol¬ 
son have noted that “a combination of expert estima¬ 
tion based on prior experience, experiments, and mod¬ 
els and simulations at virtually every mapping” is used 
in practice [206]. In many cases, these three modes of 
using experimental data are used concurrently and it- 
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Fig. 11: Schematic of information flow between experi¬ 
ments, conceptual models, and computer models. After 
Ref. [204]. 


eratively to improve the accuracy of the model and to 
bound our certainty in its predictions. 

4.1 Calibration, verification and validation 

Different communities use verification and validation 
(V&V) to accomplish different things [207,208]. While 
software developers use V&V to ensure the reliability 
and robustness of their code, the modeling community 
uses the process to quantify and improve the agree¬ 
ment between model output and experimental measure¬ 
ments (often with the additional hope that improving 
the model to address questions that are measurable ex¬ 
perimentally will help to make improvements in predic¬ 
tions to questions that are more difficult to assess). 

In building new models and simulations, we start 
first with conceptual models describing the real world 
(i.e. the experiments in the mind of the modeler). These 
conceptual models include the relevant partial differen¬ 
tial equations, initial conditions, and boundary condi¬ 
tions. The conceptual model is implemented as a com¬ 
puter model to solve these complex systems efficiently. 
Verification is the process of determining that this com¬ 
puter model is an accurate solution to the conceptual 
model. However, this process is disconnected from the 
real world, so experiments provide no role in the verifi¬ 
cation process. Observational experiments may precede 
the conceptual model by providing mechanism discov¬ 
ery (e.g. in in situ TEM) or be used to confirm the 
conceptual model is a good description of reality, but 
they do not assess whether a numerical implementation 
of a mathematical model is correct. 

The computer implementation may have adjustable 
parameters (either numerical or physical) and calibra¬ 


tion experiments may be used to set these. Only after 
a simulation tool is verified and, if needed, calibrated, 
is it validated. Validation is the process that quantifies 
how accurately the computer model represents the real 
world. Rather than asking if we solve the equations cor¬ 
rectly, it asks if we’re solving the correct equations. Val¬ 
idation is typically done through a series of high quality 
experiments under the same conditions (when possible) 
imposed by the model and it should use different data 
than any calibration data. It is often possible to break 
models up hierarchically and to validate aspects of the 
model [209]. 

The validation requirements for lower level aspects 
(e.g. nanoscale) are generally more stringent than for 
the system level aspects (e.g. continuum scale) because 
errors propagate up the hierarchy. In multiscale models, 
this often means very accurate measurements must be 
taken at very small length and time scales. Excellent 
validation has been performed using nanoindentation 
experiments, for example [210,211]. Unfortunately, val¬ 
idation at these short scales is not always easy. Available 
small length-scale models (e.g. MD) tend to also oper¬ 
ate over very short time periods and available longer 
length-scale models (e.g. CPFEM) tend toward coarser 
time scales as well. Conversely, experiments that have 
very high spatial resolution (e.g. APT) tend not to be 
able to measure fast changes and current techniques 
to observe faster kinetic processes (e.g. nanoindenta¬ 
tion) tend to be limited towards longer length scales. 
Furthermore, whereas MD may describe behavior in 
the picosecond to nanosecond regime, our time-resolved 
measurements tend to be limited to the microsecond to 
millisecond regime [212]. 

4.2 Model input 

The previous section touched on using calibration ex¬ 
periments to set adjustable parameters in a model. Ex¬ 
perimental data may enter models for another reason. 
While multiscale models can often make predictions 
from first principles, predictions can become more pre¬ 
cise or the calculations can be simpler if experimental 
input is used for scaling or to obtain material features 
that models don’t predict well. Atomic scale models de¬ 
pend upon interatomic potentials, and these may be ob¬ 
tained from experimentally-determined lattice param¬ 
eters, elastic constants, vacancy formation energy, and 
bonding energies. Potential fitting has been built into 
codes such as GULP [213] or LAMMPS [214], for exam¬ 
ple. At the mesoscale, fatigue and fracture are particu¬ 
larly sensitive to material defects. So, in fatigue mod¬ 
els [215,216], for example, crack and pore size parame¬ 
ters have been set from experimental values. 
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Regarding the use of experimentally-obtained mi- 
crostructural images as model input, there have been 
examples in the literature of using 2-D images for ei¬ 
ther 2-D mechanical property models [217] or to sim¬ 
ulate a 3-D microstructure that matches the statisti¬ 
cal properties of experimental datasets [218,219]. As 
discussed in more details in Section 2.3, full 3-D mi- 
crostructures are now being used in models. Destruc¬ 
tive analysis using serial sectioning techniques has al¬ 
lowed for the more accurate and quicker prediction of 
elastic moduli [220], the onset of plasticity, and frac¬ 
ture [221]. Non-destructive techniques, such as X-ray 
tomography have allowed identical microstructures to 
be used in both modeling and experimental determina¬ 
tion of mechanical properties [222-224]. 

Combining 3-D and nanometer-scale analysis shows 
promise and there has been recent work that uses APT 
data in MD [225,226] and DDD [227]. The APT’s im¬ 
perfect detection efficiency (modern instruments detect 
roughly half of all atoms that are field-evaporated) and 
imperfect spatial resolution have motivated a hybrid 
approach that combines atom positions determined ex¬ 
perimentally by APT with Monte Carlo simulations to 
create input suitable for density functional theory cal¬ 
culations [228,229]. 

4.3 Experimental discovery and data mining 

Because sensors, automated data acquisition systems, 
and data storage have become faster and less expen¬ 
sive, modern experimental methods generate a huge 
amount of data. Experiments have led to new discover¬ 
ies directly where new mechanisms and behavior are 
observed by a researcher who develops a conceptual 
model. But experiments now also serve as the basis for 
discovery through automated statistical analysis. The 
two challenges of data mining are first generating the 
data and then analyzing it efficiently. 

Experimental data collection: One method to gen¬ 
erate a huge quantity of data is through a combinato¬ 
rial study. Experiments can systematically vary param¬ 
eters. Even a single parameter, varied manually, can 
lead to deep insights. By varying micropillar dimen¬ 
sions cut by FIB milling, for instance, the size limits 
for plastic processes can be found [230]. But when high- 
throughput characterization techniques are used (that 
can often be conducted in parallel and/or be scripted 
to run in a serial fashion without interrupting the re¬ 
searcher), an even larger swath of variable space is ex¬ 
plored [231]. 

Data reduction and analysis: Sifting through this 
vast amount data poses a second challenge. One com¬ 
mon data reduction technique is Principal Component 
Analysis (PCA) [231,232]. Multivariate data often has 


many correlated variables that are responsible for an 
observation. PCA transforms those variables to a new 
coordinate system where the new variables are linearly 
uncorrelated and sorted from greatest-to-least variance. 

5 Open scientific questions 

The study organized by The Minerals, Metals & Ma¬ 
terials Society (TMS) entitled Modeling Across Scales: 
A Roadmapping Study for Connecting Materials Models 
and Simulations Across Length and Time Scales [233] 
identifies 30 critical gaps and limitations for materials 
modeling across multiple length and time scales and in¬ 
clude the inefficiencies in application programming in¬ 
terfaces, limitations in CPFEM analysis, the need for 
multiscale experiments to calibrate and validate multi¬ 
scale models, or the coupling of models and experiments 
for performance predictions of rare events. 

As discussed in the above sections, open questions 
lie in the different complementary domains detailed in 
this manuscript: characterizing initial microstructures 
and their temporal evolution experimentally, modeling 
coupled physical processes over time, linking different 
length and time scales, and managing large volumes 
of data. Without being exhaustive, a few examples of 
major challenges are raised. 

5.1 Gaps and limitations with current experimental in¬ 
strumentation and tools 

Measurements of stress fields within a hetero¬ 
geneous medium: Experimental visualizations of the 
displacement and rotation fields are now quite com¬ 
mon and provide important microstructural informa¬ 
tion. However, they are still limited in the sense that 
deformation is not a state variable from a theoretical 
point of view and the deconvolution between the elas¬ 
tic and plastic components of the strain field is difficult. 
An associated challenge is the measurement of the local 
stress fields with high spatial and temporal resolution. 
These measurements have to be done under mechani¬ 
cal loading, i.e. during in situ experiments. This is of 
particular importance for the stress concentration de¬ 
termination close to phase- or grain- boundaries for ex¬ 
ample; a research topic on experimental explorations 
closely linked to damage and fatigue behavior, for in¬ 
stance. 

One of the main difficulties is to combine a fine spa¬ 
tial resolution (on the order of hundreds of nanome¬ 
ters or below) with fine time resolution measurements. 
Presently, conventional XRD allows an acceptable ac¬ 
curacy for the estimation of the stresses through the 
quantification of the atomic plane spacing (see paper 
by Chang et al. published in this special journal sec¬ 
tion), but the collection time is still prohibitive. On the 
hand, synchrotron XRD, which allows 3-D local inves- 



16 


Remi Dingreville et al. 


tigation [234], is not yet accurate enough to allow such 
a measure. 

Closely related to these experimental roadblocks, 
the automation of data acquisition and analasis has to 
be developed to analyze many diffraction patterns in a 
relatively short time. Such improvement is critical be¬ 
cause the required accuracy is very high and the signal 
profiles and maxima are not easy to define without any 
ambiguity. A huge volume of data will be generated 
from which we must extract the most important infor¬ 
mation (see Section 5.3). The improvement of such mea¬ 
surements will not only offer more experimental data on 
complex materials systems more easily and contribute 
to the V&V of models, but it will also drive the devel¬ 
opment of more complex anisotropic models. 

Improvement of the accuracy and sensitivity of 
non-destructive 3-D measurements: 2-D investi¬ 
gations and destructive 3-D analysis examining 2-D sec¬ 
tions may introduce a bias due to free surfaces where 
mechanical and environmental states are quite different 
from the bulk. These effects can generally be accounted 
for in simulations, but the actual mechanisms may still 
be misunderstood. While computational power and al¬ 
gorithms are allowing 3-D simulations that yield greater 
insights, non-destructive 3-D characterizations must be 
made available with a similar accuracy and sensitiv¬ 
ity. During the past ten years, many new developments 
have been made through synchrotron radiation. Some 
recent examples in the literature demonstrate local mis- 
orientation measurement [131] and the localization of 
nano-cracks by X-ray tomography [56]. 

Improvements can also be made in the speed of some 
experimental measurements, which tend to be much 
slower than desired strain rates. Ideally, the measure¬ 
ments would be coupled with in situ deformation and 
could be done on the fly as the deformation occurs. This 
implies the development of mechanical experiments on 
small samples deformed with devices compatible with 
the diffraction and other investigation tools. The lim¬ 
itation will probably be defined by the representative 
volume which allows a continuous mechanical descrip¬ 
tion (see Section 5.2). 

Combination of simultaneous measurements at 
the same location: An additional challenge now is 
to compile different types of 3-D measurements of the 
same sample (ideally, simultaneously) to determine lo¬ 
cal microstructure-properties relations. There are now 
efforts to correlate APT data with electron tomography 
and FIB measurements to improve reconstructions, for 
example [235]. 

A greater ability to combine information from mul¬ 
tiple analysis techniques and a greater ability to per¬ 


form the techniques simultaneously would be powerful 
additions to our toolbox. Ideally, one would be able to 
measure not only chemical segregation, but also corre¬ 
late that with the knowledge of local strains (or elastic 
stresses), particularly close to phase or grain boundaries 
where segregation can play an important role for me¬ 
chanical localization. Similarly, our models may benefit 
if grain character could be correlated with local yield 
strength. 

5.2 Gaps and limitations with current models 

Scale multiplicity due to microstructural evo¬ 
lution and materials response: The core principle 
that materials properties are determined by composi¬ 
tion and structure requires that predictive models are 
able to treat multiple scales within the same paradigm 
consistently. While many macroscopic behaviors result 
from the average of local variables, such as the effect of 
grain size or phase fraction on the stress-strain curve, 
some are directly linked to specific localizations that 
result from microstructure distribution and the chance 
of deformation, impurity diffusion, phase changes, etc. 
This dependence on localized quantities is the norm 
for damage, fracture, and fatigue. Because these may 
arise in a large structure, one has to combine a com¬ 
plete macroscopic description that accounts for bound¬ 
ary conditions with an accurate quantification that al¬ 
lows local calculations. For a nanosized sample, the 
characteristic volume of the microstructure is of the 
same order of the representation volume, but for meso- 
or macro- scopic samples, a complete description of the 
microstructure at a small scale is impossible. How can 
we overcome this problem? Periodic structure assump¬ 
tions and/or homogenization are possible, but are not 
accurate for local phenomena such as damage, cracks, 
precipitates, or recrystallization nuclei. There is a real 
theoretical challenge to represent the whole microstruc¬ 
ture with its average heterogeneity. Current models em¬ 
ploy discretization methods on a specific grid size and 
time marching strategy. The major challenge associ¬ 
ated with dealing with disparate scales within the same 
framework resides in the fact that different scales re¬ 
quire different grid resolution and time steps, making 
the scale bridging a complicated endeavor. A spatial 
example includes first- and second- generation twins in 
alloys [107, 236, 237] while a temporal example is the 
phenomena associated irradiation with ultrafast defect 
generation (displacement cascade) versus the slow de¬ 
fect migration and interaction [238]. Such considera¬ 
tions will undoubtedly need to be developed in con¬ 
junction with the development of new instrument while 
expanding their accuracy and sensitivity across scales. 
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Quantification and propagation of uncertainties: 

Many fundamental processes and mechanisms, when 
viewed at the appropriate scale, behave stochastically. 
Atomic diffusion is one example and yet the collective 
behavior of large numbers of diffusing atoms can, on 
average, be quantified deterministically. Materials com¬ 
position too can, depending on the scale, be consid¬ 
ered highly heterogeneous. Uncertainty quantification 
and its propagation across scales (time and space) are 
both very challenging. Indeed, uncertainty at different 
scales needs to be considered differently. The increas¬ 
ing cross-pollinations between the modeling community 
and the experimental community will also need to bring 
in the statistics and probabilistic communities in order 
to address the variability and predictability of mate¬ 
rials response over multiple scales. Integrated rigorous 
statistical paradigms will ultimately permit the identifi¬ 
cation of critical microstructural features that are vital 
to property evolution, especially at the system level, 
leading to materials design solutions. 

Inefficiencies of current approaches for multi¬ 
scale modeling: A wide range of multiscale model¬ 
ing approaches is limited by the increasing computa¬ 
tional expense and cumbersome algorithms associated 
with dealing with multiple scales and/or large sets of 
experimental data that need to be incorporated into 
the models. This include pre-conditioners, i.e. reliable 
mathematical techniques conditioning the input data 
and problem into a more suitable form for numerical 
simulations or the improvement of reduced order mod¬ 
eling methods (i.e. reducing the degrees of freedom) 
enabling limitation of resolution losses. An additional 
motivator to address these inefficiencies is the promise 
that early model results may help to design a set of 
experiments that explores critical parts of the parame¬ 
ter space rapidly and inexpensively, leading to further 
model improvement. 

5.3 Gaps and limitations with data management 

Storage and sharing of “big data”: The volume of 
data (both from numerical models and experimental 
measurements) that will be acquired within this emerg¬ 
ing field of Materials Science is undoubtedly huge. The 
heterogeneous, complex and convoluted nature of the 
data collected will need to be distributed across many 
communities (see Section 1). Resolving the challenges 
associated with the data management will require the 
establishment of common standards and protocols to 
vet, compare and use such data reliably. Such mecha¬ 
nisms need to be put in place for the community to not 
only allow to compare, verify, and validate data used 
in models, but also to provide means to make decisions 
based on greater amounts of information entering the 


model (with no need for recalculating/charactering). 
Kalidindi and De Graef [144] provide a more detailed 
review on important aspects of materials data manage¬ 
ment, such as storage hardware, archiving strategies, 
and data access strategies. Marsden et al. [239] high¬ 
light the breadth of data that must be stored for mul¬ 
tiscale materials modeling and the challenge of making 
this data interoperable with many computational tools. 

Speed up engineering transfers of fundamental 
approaches: Even though Materials Science is con¬ 
tributing to knowledge development, its final goal is 
to provide technological innovations and practical tools 
and knowledge relevant to the industry. One of the main 
practical advantages of multiscale approaches is orient¬ 
ing new material developments through an inverse ap¬ 
proach. That is, knowing the exact distribution of final 
properties expected in a part, to be able to design the 
optimum material, combining graded properties and 
multi-material solutions. Such an approach gives rise 
to guidelines for chemical composition and processes 
for new alloys and materials. However, it is untenable 
to implement complex multiscale, multi-physics codes 
in an industrial context. It is then of importance to 
develop simplified simulations integrating most of the 
specificity of the fundamental codes. A complementary 
method is to design coupled approaches for which dif¬ 
ferent codes interact at different scales. For instance, 
a 2-D finite element code with macroscopic behavior 
laws may identify the most critical part during forming 
of a complex part, generates the boundary conditions 
to a specific volume, and a more complete 3-D physics- 
based model can probe the part. Some of these methods 
are already in place, but some work has to be done to 
improve interactions between the different scale models 
and to speed up calculations. Finally, databases have to 
be created and qualification of already existing materi¬ 
als is needed in order to reduce the domain that models 
would explore. 

Inefficiencies of programming interfaces: Porta¬ 
bility between different codes, algorithms, and instru¬ 
ment interfaces are not always available and/or are in¬ 
efficient because of computational complexity to map 
data into a usable form. For example, phase field models 
coupled with thermodynamic databases for the multi- 
component systems lack the ability to interface effi¬ 
ciently and information transfer from one to the other 
is often cumbersome and amateur. The challenge asso¬ 
ciated with this interfacing problem is to devise bidi¬ 
rectional methods and protocols to efficiently commu¬ 
nicate and pass data. This statement is true not only 
for interfaces between two models but also for models 
directly incorporated with experimental setups [203]. 
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6 Concluding remarks 

Multiscale modeling approaches have proven to be pow¬ 
erful methods to predict macroscopic mechanical re¬ 
sponses for a variety of loading conditions. They rely 
on experiments for validation, model input, and mech¬ 
anism discovery. As these approaches are extended to 
new problems—new material systems, environments, 
and/or forms of mechanical response—new models, al¬ 
gorithms, and experimental techniques are needed and 
will also need to be integrated within the same frame¬ 
work in order to meet the increasing demand for better 
insights and predictability of the performance of such 
materials systems. 

Recent developments in experimental instrumenta¬ 
tion have helped build better multiscale models. Break¬ 
throughs in materials research instrumentation now en¬ 
able us to link atomistic, nanoscale, microscale, and 
mesoscale phenomena across time scales. Emerging re¬ 
search opportunities for experimental improvements in¬ 
clude increasing resolution and fidelity, decreasing mea¬ 
surement time, automation of analyses, and better man¬ 
agement of the large amounts of data. Such improve¬ 
ment will undoubtedly lead to a much deeper under¬ 
standing of the local phenomena associated with hetero¬ 
geneities and interfaces. Finally, 3-D investigations are 
now largely able to quantify microstructural evolutions 
in bulk and not only on the sample surface with the 
influence of free surfaces, or on a section made though 
destructive preparation. 

Along with these experimental advances, numerical 
and computational improvements enable the study of 
complex materials systems across multiple length and 
time scales. Even though the available experimental and 
computational tools provide incomplete and uncertain 
information on phenomena of interest, improvements 
in the affordability of computational power and storage 
capabilities allow the simulation of large representative 
volumes capturing more and more microstructural fea¬ 
tures. New algorithm developments speed up the calcu¬ 
lation and lead towards time scales of simulations that 
are compatible with experiments. One of the main chal¬ 
lenges associated with these models is related to con¬ 
necting information across scales to reduce the amount 
and complexity of the information at higher scales: we 
have no a priori knowledge of which physics and vari¬ 
ables at the lower length scales are “worth” upscaling. 
As such, emerging research opportunities for modeling 
improvements include novel approaches to identify this 
information. 

Finally, both of these improvements are accompa¬ 
nied by an increasing focus on materials data analyt¬ 
ics, which focuses on data-driven approaches in Mate¬ 
rials Science to extract and curate materials knowledge 


from available data sets (from numerical and experi¬ 
mental studies). Opportunities in this held of research 
include the standardization and automation of work- 
flows to handle “big data” and develop protocols and 
strategies to take advantage of emerging high perfor¬ 
mance data science toolkits and cyber-infrastructure 

Despite the impressive progress made during the 
last decade, the number of materials features describing 
the links between processing, structure, and properties 
is too large to allow systematic study using established 
scientific approaches with modest resources. The Mate¬ 
rials community is still facing huge challenges as legacy 
data, numerical models, and multi-physics knowledge 
are fragmented and compartmentalized across so many 
disciplines. Research at the crossroads between com¬ 
putational modeling and experimental characterization 
will continue to have a broad impact on a wide range of 
science and engineering problems in Materials Science, 
but only through an efficient and interactive interdis¬ 
ciplinary collaboration. Given this trend, the develop¬ 
ment of modern data science and cyber-infrastructure 
is critical to mediate and accelerate such collaborations. 
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